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1.  Introduction 


Advanrrs  in  molecular  beam  epitaxy  and  other  crystal  growth  t('ch- 
niques  have  allowed  detailed  studies  of  semiconductor  heterostructures 
[1].*  In  the  last  few  years,  the  double-barrier  resonant  tunneling  struc¬ 
ture  has  been  studied  by  a  number  of  groups.  This  structure  consists 
of  a  thin  layer  of  CaAs  (about  50  A)  surrounded  on  both  sides  by  a 
layer  of  .MGaAs  (of  similar  thickness),  all  embedded  in  a  Ca.As  ( rystal. 
When  electrically  contacted  and  biased,  this  structure  disi)lays  peaks 
in  its  current-voltage  characteristic  [2-13].  From  the  scattering  point  of 
view,  these  peaks  are  attributed  to  peaks  in  the  transmission  coefTicicMit 
[14].  In  terms  of  the  electronic  structure,  the  double-barrin-  potential 
confines  electronic  states  of  certain  energies,  creating  re.sonances  that 
manifest  themselves  as  peaks  in  the  one-dimensional  local  density  of 
states  between  the  barriers  [15-18].  For  a  symmetric  stiuictiire  in  zero 
applied  electric  field,  the  states  can  be  labelled  by  even  or  odd  parity, 
o  =  c  or  o.  a  wave  vector  transverse  to  the  growth  direction,  k,  and  a 
wave  vector  along  the  growth  direction,  q.  associated  with  the  solution 
of  a  one-dinuMisional  effective  ma.ss  Schrodinger  equation.  The  doubly 
degenerate  electron  energy  eigenvalues  are  given  by 

p:  =  h\]^ +  q^)/2wu  (1) 

where  rrti  is  the  conduction  band  effective  mass  at  the  F-point  (for 
CaAs).  For  a  given  value  of  k,  there  are  resonant  values  of  q  at  which 
the  amplitude  of  the  wave  function  inside  the  cjuantum  well  region  is 
large  compared  with  the  wave  function  amplitude  outside  this  region. 
For  wave  vector  vahvs  q  that  arc  not  close  to  the  resonant  values, 
the  amplitude  of  the  wave  function  inside  the  quantum  well  region 
is  suppressed  relative  to  its  outside  value.  Outside  the  quantum  well 
region,  the  wave  function  amplitude  is  roughly  the  same  for  energies 
on-  and  off-resonance. 

The  above  description  of  electronic  states  neglects  certain  aspects 
of  the  problem,  such  as  external  electric  fields,  scattering  du('  to  ])ho- 
nons  [19],  im|mrities.  and  electron/electron  interactions.  The  ('fft'cls  of 
electron/electron  int«'ractions  (space-charge  buildup)  on  current -volt age 
characteristics  have  been  discussed  [20]  and  treated  in  rletail  within  a 
self-consistent  scheme  [21].  Inelastic  scattering,  and  its  consequences 
for  resonant  tunneling,  has  also  been  discussed  by  sev('ral  authors  [18,22]. 
However,  elastic  scattering  due  to  imiuuiti('s  has  not  bcxMi  tr('at<'d  the¬ 
oretically.  des]-)itethe  fact  that  at  licpiid  nitrogen  temiunat  ures,  where 
measurements  on  double-barrier  structures  are  freriuently  made,  (dastic 
'References  are  listed  at.  the  end  of  the  main  body  of  te.xt. 
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scattering  caused  by  charged  impurities  is  the  predominant  scattering 
mechanism  [23].  While  it  is  true  that  elastic  scattering  events  are  not 
expected  to  have  an  effect  on  resonances  that  is  comparable  to  inelastic 
events,  recent  experiments  [13]  on  intentionally  doped  double-barrier 
structures  show  that  their  effect  is  nontrivial  and  needs  to  be  better 
understood.  In  this  paper  we  examine,  in  a  qualitative  way,  the  man¬ 
ner  in  which  scattering  rates  from  a  single  charged  impurity  depend 
upon  the  electron’s  initial  state  and  the  impurity’s  position.  Next,  as¬ 
suming  there  is  a  distribution  of  randomly  positioned  impurities,  we 
add  the  probability  for  scattering  due  to  each  impurity,  resulting  in  a 
total  elastic  scattering  rate.  We  find  that  for  a  finite  system  size  and 
certain  values  of  the  parameters,  such  as  impurity  potential  range,  the 
scattering  of  resonant  electrons  (in  the  limit  of  low  electric  fields)  can  be 
dominated  by  impurities  in  the  quantum  well  region.  This  is  surprising 
because  in  circumstances  where  sharp  resonances  are  absent,  contribu¬ 
tions  from  the  quantum  well  region  arc  the  same  order  of  magnitude  as 
surface  effects,  which  can  be  ignored. 

In  section  2  we  establish  our  notation  and  briefly  summarize  the 
form  of  the  electron  eigenstates  for  a  double-barrier  potential.  In  sec¬ 
tion  3  we  use  these  states  to  compute  the  scattering  rate  due  to  one 
impurity,  and  in  section  t  we  consider  the  effects  of  a  random  distribu¬ 
tion  of  impurities. 

2.  Eigenstates 

Within  a  one-band  effective  mass  theory,  an  electron  state  is  represented 
by  a  product  wave  function  [24] 

'l'c,.k.,(r)  =  i/’a,k,,(r)nc(r)  =  (r|o.k.(?),  (2) 

where  ^tc(r}  is  the  conduction  band  Bloch  function  at  the  F-point  (and 
is  taken  to  be  the  same  function  for  Ca.'\s  and  .AlCa.As),  and  ra.k.,]*") 
are  the  envelope  functions.  For  a  double-barrier  potential 

‘'-'"lo  otherwise, 
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the  envelope  functions,  normalized*  to  unity  over  a  volume  .  wf'rc 
previously  calculated  [17]  and  can  be  written 

=  ^exp(?;k  •  r)4>„.k,,(-).  (1) 

where  k  =  and  q  labels  the  parity  of  the  state,  q  =  f(o)  for 

even  (odd)  parity.  The  ^-dependent  functions  are  given  by 

where  6a,e  and  6^,0  are  Kronecker  delta  functions,  and  the  auxiliary 
functions  d^in),  c„(?/),  fa{u),  5fo(n),  and  F^,k(^0  are  defined  in  ap¬ 
pendix  A.  The  uinction  w{u)  is  given  by 

(6) 

and  carries  a  dependence  on  transverse  wave  vector  k  through 

7='^-(l-/?)lklV.  (7) 

Here,  7  is  the  dimensionless  barrier  height  given  by  7  =  2mi\oa^//j^, 
and  /?  is  the  ratio  of  the  effective  mas.ses,  /?  =  »77i/m2,  where  mi  is  the 
mass  in  the  quantum  well  region  and  outside  the  barriers,  and  mj  is  the 
mass  in  the  barrier  region.  Because  the  effective  mass  in  the  quantum 
well  region  is  different  from  that  in  the  barrier  regions,  the  effective 
barrier  height,  7,  is  different  for  each  k. 

The  functions  and  are  continuous  at 

the  interfaces  z  =  d:a  and  z  =  dz{l  -{-S)a.  Periodic  boundary  conditions 
in  a  box  of  area  Ij  x  L  have  been  used  in  the  x  and  y  directions, 
so  the  allowed  values  of  k  =  2:r(7?j.,  n^,  0)/L,  where  iij.  and  Vy  can 
take  any  integer  values.  It  is  interesting  to  note  that  at  sufficiently 
large  |k|,  the  effective  barrier  height  may  actually  be  negative.  The 
allowed  values  of  the  ^-component  wave  vector,  q,  are  given  by  the 
roots  of  a  transcendental  equation  whose  prc'cise  form  depends  on  the 

'In  normalizing  the  wave  function,  we  have  dropped  terms  of  0(n/f.)  and 
0{a6/L).  In  tlie  quantum  well  limit,  7  — *  oc'  or  6  --  oc',  the  function  /o(7<)) 
has  zeros  at  re.sonant  values  of  q,  and  the  term.s  0(a/L)  and  0{aS/L)  cannot  he 
dropped.  However,  for  any  finite  value  of  7,  there  is  always  a  large  enough  L  so 
that  dropping  these  terms  is  justified.  Physically,  this  means  we  are  neglecting 
surface  effects,  which  can  l)e  important  when  the  thickness  of  the  crystal  hecome.s 
sufficiently  small,  or  when  W'e  approach  the  7  —  oc  or  ^  *  00  limit,  with  L  fixed, 
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k(<?«)  f  +  sin(<72)4,o  |c|  <  a 

'~f=^  I  dj^q(i)ey^{w{qa)zla)  +  ea{qa)  e\p{-w{qa)z  j a)  a  <  z  <  {\  -\-  8)a  (5) 

i  fo,iqa)(;xp(iqz)  +  gc{qa)exp{-iqz)  (1  +  <  ^  <  f/2 


boundary  conditions  at  ^  =  ±Z//2.  However,  in  the  limit  aj L  — ►  0,  the 
density  of  roots  approaches  L/2;r,  and  the  solution  of  this  equation  is 
not  necessary.  The  energy  eigenvalues  associated  with  the  state  given 
by  equation  (2),  measured  from  the  bottom  of  the  conduction  band, 
are  given  by  equation  (1). 

The  envelope  functions  given  in  equation  (4)  contain  the  resonance 
effects  due  to  the  presence  of  the  double-barrier  potential  through  the 
functions  Fgy{u)  and  Foy,{u).  For  a  given  k,  these  functions  have  peaks 
at  u  =  qa,  for  wave  vector  values  q  =  ge,i>  9e,2)  9e,3,  •  •  and  qo,\,  qo.2' 
qo,3^  .  ■ respectively.  The  peaks  in  these  functions  correspond  to  an 
even  and  an  odd  series  of  resonances.  These  functions  are  intimately 
related  to  the  local  density  of  states,  and  have  been  discussed  in  detail 
elsewhere  [17].  An  electron  with  a  given  transverse  wave  vector,  k, 
and  q  close  to  a  resonant  value  has  a  wave  function  with  a  large  am¬ 
plitude  in  the  quantum  well  region  and,  correspondingly,  a  decreased 
amplitude  outside  the  barriers,  so  that  normalization  over  the  interval 
(  — L/2, -t-L/2)  is  maintained.^  An  electron  whose  wave  vector  is  not 
near  resonance  has  a  small  wave  function  amplitude  in  the  quantum 
well  and  barrier  regions  and,  correspondingly,  a  larger  amplitude  out¬ 
side  these  regions.  Figure  1  shows  a  plot  of  the  square  of  the  even 
envelope  function,  |V’e,k.7(r)P,  versus  z/a,  for  on-resonance  and  off- 
resonance  electron  states.  The  state  that  has  a  wave  vector  in  the 
resonance  region,  qa  =  0.80,  has  a  large  v.ave  function  amplitude  in  the 
well,  compared  to  a  state  whose  wave  vector  lies  outside  the  resonance 
range,  qa  =  1.5. 

When  the  resonances  are  quite  sharp  it  is  appropriate  to  speak  about 
resonant  subbands,  analogous  to  the  two-dimensional  subbands  in  a 
quantum  well  structure.  A  resonant  subband  is  defined  by  a  re.sonant 
value  of  q  (equal  to  one  of  the  q^,:,  qe.2,  ?e,3,  •  ■  or  (?„,],  90,2,  90.3.  •  •  •) 
and  a  continuous  value  (in  the  limit  L  — >  00)  of  transverse  wave  vector 
k.  These  resonant  subbands  differ  from  those  in  a  quantum  well  in 
an  essential  way.  In  a  quantum  well  the  subbands  arise  from  truly 
localized  states  (for  energies  below  the  top  of  the  barriers),  while  in 
a  double-barrier  structure  the  resonant  subbands  arise  from  e.xtended 
states. 

hSee  footnote  on  page  7. 
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Figure  1.  Square  of  even  function,  L^|T/’e,k,(}(r)!^,  is  i)lottecl 
versus  spatial  position  zja^  for  “on-resonance” 

{qa  =  q^  ia  =  0.80),  and  “off-resonance”  {qa  =  1.5),  electrons  for 
case  k  =  0,  ^  =  1 .0,  7  =  2.0,  and  13  =  0.73. 

3.  Scattering  Rate  Due  to  One  Impurity 

The  elastic  scattering  rate  out  of  a  conduction  band  state  (ft,k,('/),  to 
any  other  conduction  band  state,  can  be  found  from  Fermi’s  golden  rule 

=  x  ^  |(a',k',<7'K/|a,k,,/)p<^(Ek,,-Ekv)’  («) 

”  a'.k'.?' 

where  E(r)  is  the  impurity  potential.^  For  many  impurity  potentials 
the  matri.x  element  in  equation  (8)  can  be  approximated  as  a  matrix 
element  between  effective  mass  wave  functions  (see  app  B): 

(a',k',^'Ky|o,k,9)=  /  dV:<,kv(r)E(r)7’„,k,,(r).  (9) 

Calculation  of  this  matrix  element  can  be  further  simplified  by  substi¬ 
tuting  into  equation  (9)  the  explicit  form  of  the  wavefunction,  given  in 
equation  (4),  and  using  the  Fourier  expansion  of  the  potential, 

b'(r)  =  X!E(q)exp(?q  •  (r  -  R,,)).  (10) 

q 

*  VVe  arc  not  ron.siOering  potentials  U{r),  which  couple  to  spin  degrees  of  freedom; 
hence  we  have  omitted  the  spin  labels  on  the  electron  states.  The  cpiantily  ^ 
is  the  scattering  rate  out  of  a  state  labelled  by  o,k,7  with  either  spin  up  or  down. 
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where  R^,  is  the  impurity  position.^  The  x  and  y  spatial  integrations  can 
then  be  performed,  leading  to  a  two-dimensional  delta  function  in  wave 
vector  variables.  The  sum  over  the  x-  and  y-wave  vector  components 
can  then  be  done,  resulting  in  the  following  expression: 

(o',k',9'K/|cv,k,y)  =exp(-*{k'-k)  -RJ 

J-L/2 

•^f'/(k'  -  V->r  p^z)c\p{ip,{z  -  Zo)),  (11) 

Pz 

where  Zg  is  the  ^-component  of  the  impurity  position  vector  and  z  is 
a  unit  vector  along  the  growth  direction.  The  sum  on  p^  is  over  the 
^-component  of  the  wave  vectors  appearing  in  equation  (10),  which 
were  u.sed  to  expand  the  impurity  potential.  The  x-  and  y-components 
of  the  impurity  position  vector  in  equation  (11)  enter  only  as  a  phase. 
Because  the  system  is  translationally  invariant  in  the  x  and  y  directions, 
there  is  no  loss  in  generality  if  we  choose  the  x  and  y  components  of 
the  impurity  position  vector  to  be  zero. 

Using  equation  (II)  we  can  immediately  make  a  rough  estimate  of 
the  ratio  of  probability  for  scattering  of  re.sonant  electrons  within  the 
first  subband,  |e,k,  ye.i)  '■o  scattering  out  of  the  first 

subband  into  the  second  subband,  |e,k,ye.i}  — >  1®,  k',y„,i),  by  noting 
that  the  function  proportional  to  This  ratio  is 

proportional  to  the  factor^ 

^<f,k'(ge.lQ) 

^o,k'(9o.i^) 

and  is  independent  of  the  detailed  form  of  the  impurity  potential.  Here 
qe,\  and  y'  ,  are  the  resonant  wavevectors  corresponding  to  transverse 
wave  vectors  k,  and  k',  respectively.il  For  the  structure  used  by  Sollner 
et  al  [3],  this  ratio  is  =^35,  for  k  =  k'  =  0. 

^We  assume  that  the  impurity  is  located  far  from  the  walls  of  the  box  compared 
to  any  length  scale  of  interest. 

'This  result  is  true  only  for  impurities  located  in  the  quantum  well  region.  Equa¬ 
tion  (5)  may  suggest  that  at  re.sonance  the  wave  function  is  large  beyond  the  barri¬ 
ers,  r  >  ( 1  ^)a,  because  Jj(?<i)  is  large.  This  is  not  the  case  because  the  factor 

/a('7<i),  which  multiplies  j{(<?a)  in  the  wave  function,  becomes  small  at  resonance, 
leading  to  a  wave  function  that  is  not  large  outside  the  quantum  well  region. 

ill  he  re.sonant  wave  vector  values  q  =  q^  i,  7^  2,  q  =  7,  3,  . .  .,  and  q^j,  qo,-2,  qo,3< 
.  .  .,  depend  mildly  on  the  corresponding  value  of  k.  This  is  a  consequence  of  the 
different  effective  masses  in  the  quantum  well  and  barrier  regions  (see  eq.  (6)). 
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In  order  lo  get  a  elearer  picture  of  (‘!a.stic  scattering,  we  model  a 
charged  impurity  by  a  spherically  syiimietric  scre<-ned  Conloml)  i)oten 
tial  [-ire-i!)]: 

K  I  r -  R„  I 

where  1/A  is  the  range  of  the  pot('ntial,  Ze  is  the  iir.pnrity  charge,  R,, 
is  (he  impurity  position,  and  k  is  the  static  dielect ic  constant.  I'or 
simplicity,  we  are  neglecting  the  effects  of  anisotropic  scrcs'iiing.'*  for 
this  impurity  potential,  the  sum  on  in  equation  (11)  can  b<'  changed 
to  an  integral  and  evaluated,  resulting  in  the  form 

Jtt  Z(  ^  1 

(o',k',»/'|f ’lo.k.f/)  =  ^ - /  (I'M 

IZ  K  S  J-I.Ji 

where  the  effective  scr<<'ning  w'ave  v<'ctor,  s,  dc'pends  on  tlu'  transverse 
moiiKMitum  through  tlu'  r<'lation 

s'^  =  A' +  (k  -  k')'.  (I  I) 

Using  ('(juat  ion  ( 1  d),  we  define'  tlu'  dime'usionh'ss  (piant  ity  /,,',k'.7':.>,k,v(  -.d 
by 

{<i',k',</|r|(,,k.</)  5  - - ! - ;T7J  r.',k'., (I'i) 

Aye  +  (k  -  k')V^ 


wh('r('  we  have  introduced  the  Bohr  radius  <i„  =  s  j  in  ^  Z(  ^ .  (’alcu- 
latieui  of  the  function  /„<,k',,,';,v,k,i)  tedious  but  straightforward,  and 
tli<'  results  for  an  impurity  locate'd  inside  the  (piant um  wc'll  rt'gion  and 
inside  the  iuirrier  regions  are  given  in  ajipendix  ('. 

In  the  special  ca.se  when  the  im|)urity  is  locati'd  insidi'  the  wi'll,  or 
beyond  the  barriers  (;„  >  (1  y  A)(i),  and  when  .s  is  sulliciently  huge 
i.e.,  |k  —  k'l"  is  large  (but  still  (k  —  k')-c  <C  1,  where  c  is  a  primitive' 
lattice'  ve'ctor)  the'  matrix  e'le'me'iii  can  be'  a|)pie)ximat('d  by 

,,,,,,,  ,  .  Sir  Ze‘ {/■!,'.k'(</'«)/''.,.k'//")}  ^  ,  - 

(o  .k.e/  |f  |o,k,,/)  - - (k~- k')'^ - '^^:.^k^,'(-)‘^,.,k, (lb) 

where'  W('  ha\e  assume'd  .s/>  ^  1,  droiiping  te'rms  of  oreh'r  ('xp( —s/,/2). 

I'his  re'sult  show's  that  the  matrix  e'le'ine'nt  is  large'  at  re'sonaut  \alu('s 
of  initial  and  final  e'le'ctron  wave'  ve'ctors.  wlu're'  /',,',k'('/^" )  ‘'’'d  /'l..k('/^') 

"A  corrcfl  trratiiK'iil  of  srp'ciiing  in  a  (ioiil'lr  harrirr  slrihiuri'  is  a  (litluiill 
prol'lf'Mi  in  ilsi’lf  WV  havr  rliosi-n  (<<  inoili'l  tin'  impurity  t'v  a  phvsifallv  plaiisihli 
/(TdI Ii  onliT  poti'iilial  with  a.s  few  unknown  parameters  .as  possihli'  Int roifuring 
anisotro|)if  screening  would  iimiouhtedlv  introduce  mi'ri-  unknown  paranu'tf  rs 


^  ,.ve  peaks.  As  |k  —  k'|  increases,  the  magnitude  of  the  matrix  ele¬ 
ment  decreases.  This  is  a  consequence  of  the  rapid  oscillations  in  the 
transverse  part  of  the  wave  function  that  serve  to  diminish  the  integral. 

We  now  substitute  the  explicit  form  of  the  matrix  element  given  in 
equation  (15)  into  equation  (8)  for  the  scattering  rate  and  sum  over 
all  final  states.  In  the  limit  L  — »  oo,  the  density  of  wave  vectors  ap¬ 
proaches  (Lf'lTr)^  per  unit  volume  in  (k,^)-space,  and  we  can  change 
the  sum  over  wave  vectors  to  an  integral  over  the  region  q'  >  0  and 
—  oc  <  <  +00.  In  cylindrical  coordinates,  we  do  the  q  integra¬ 

tion,  leading  to  the  following  expression  for  the  scattering  rate: 


a 


(17; 


where 


,  f-l-l 

=  I  dk'k'J^ 


dx- 


[k'^  +  q^  -  k'-^yi^ 

■  ^  +  G(a',F,7;Q,I-,f/; -x)].  (18) 

a'=e.o 

In  equation  (18)  we  have  introducecl  the  dimension le.s.s  (unction 


C;{a'.k'.q\oiAk,q\x)  =  ^2^2  q.  (k  _  kTy^a^' 

and  the  dimensionless  wave  vectors  k  =  |k|a,  k'  =  |k'|fi,  q  =  qa.  la 
the  integrand  of  equation  (18),  the  functions  G  are  -valuated  at  q'  =  q, 
where  q  =  [A  +  q^  ~  The  integration  is  over  the  dimensionless 

final  transverse  wave  vector,  k',  and  x  =  cos(d)),  where  0  is  the  angle 
between  the  initial  and  final  wave  vectors,  k  and  k'.  I'he  integration 
over  k'  has  a  finite  upper  limit,  a  consequence  of  energy  conservation. 
I'he  double  integral  in  equation  (18)  was  done  numerically,  and  the 
results  are  shown  in  figures  2  ;  '  d  3.  for  parameters  characterizing  the 
structure  used  by  Sollner  et  al  [3].  for  an  im])urity  potential  range 
1/A  =  lOn. 

'fhe  solid  curve  in  figure  2  is  the  scattering  rate  out  of  an  (wen  state, 
riu'  peak  at  qa  =  0.898  (off  the  scale)  is  due  to  the  lowest  energy  evi'ii 
resonance,  while  the  peak  near  qa  =  2.0  is  due  to  the  next  higher  en¬ 
ergy  even  resonance.  The  dashed  line  shows  the  scattering  rate  out  of 
odd  states,  and  the  peak  at  qa  =  1.66  is  due  to  the  lowest  energy  odd 
resonaiux'.  ,'\s  one  would  exp<  ct,  the  peaks  in  tlu'  even  and  odd  scatter¬ 
ing  rates  alternate  with  increasing  wave  vector.  Imr  this  structure  and 


potential  range,  the  ratio  of  the  scattering  rate  at  the  largest  ev<>n  |)eak 
to  the  largest  odd  peak  is  approximately  100.  The  dottc'd  line  in  figure 
2  shows  the  scattering  rate  out  of  even  states  in  a  hulk  semiconductor 
(where  7  =  0  and  /?  =  I)  with  the  same  range  potential  1/A  =  lOo. 
The  large  scattering  rate  near  g  =  0  is  due  to  the  Coulornhic  nature 
of  the  impurity  potential  (there  would  be  a  divergence  for  1 /A  — ►  oc). 
Electrons  which  occupy  states  near  the  lowest  energy  resonance  expe¬ 
rience  scatterings  3000  times  more  frequently  than  occurs  in  the  hulk. 
For  a  longer  range  potential,  the  peaks  in  the  scattering  rate  would  he 
less  pronounced.  However,  the  overall  scattering  rate  would  he  bigger 
because  of  increa.sed  wave  function  and  impurity  potential  overlap.  Be¬ 
cause  of  the  large  scattering  rates  for  states  near  resonance,  we  expect 
that  the  local  density  of  states,  averaged  over  the  quantum  well  region, 
can  he  significantly  modified  from  that  calculated  without  impuritic-s 
[15,17]. 
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0.0  0.5  1.0  15  2  0  2  5  3.0 
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Figure  2.  Dirnensionle.ss  .scattering  rate, 

(8/(«'^),  is  plotted  versus  qa,  for  an  impurity  located  at  center  of 
quantum  well,  Zo  =  0.  Solid  and  dashed  curves  are  scattering 
rales  out  of  even  and  odd  states,  a  =  c  and  rv  =  o.  respectively, 
for  k  =  0  and  harrier  height  V/,  =  0.23  eV,  effective  mass  in 
quantum  well  mj  =  0.067,  6  =  2.0,  P  -  0.73,  and  impurity 
potential  range  1/A  =  10a.  For  case  o  =  r,  p<'ak  at  qn  =  O.SOS  is 
off  the  scale,  reaching  a  maximum  value  of  ~  3000,  in  units  on 
graph.  Dotted  curve  is  scattering  rate  in  ahsi'uce  of  a 
douhle-hariT'r  structure,  7  =  0.0,  p  —  1.0,  for  a  pot('ntial  rang(' 
1/A  =  10a. 
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F  igure  3.  Solid  line  is  a  plot  of 
dimensionless  scattering  rate. 

( TO  1  Aq  />^ ) / (8/ia^ ) ,  versus 
impurity  position  Zo/a,  for 
resonant  state  k  =  0, 

=  0.897,  for  same  structure  as 
in  figure  2.  Dashed  line  below 
solid  curve  shows  position  of 
barrier. 
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Figure  3  shows  a  plot  of  the  scattering  rate  as  a  function  of  impurity 
position  for  states  in  the  range  of  tlie  lowest  energy  even  re.sonance. 
Impurities  near  the  center  of  the  ejuantum  well  give  a  somewhat  larger 
rate  of  scattering  than  those  near  the  edges,  because  of  the  peaked 
nature  of  the  even  wave  functions.  However,  for  this  potential  range, 
this  is  a  minor  effect.  In  the  calculations  above,  we  made  a  somewhat 
arbitrary  choice  of  the  potential  range;  1/A  =  10«.  FAir  a  shorter  range 
potential,  the  dependence  of  the  scattering  rate  on  impurity  jjosition 
would  be  stronger  than  that  shown  in  figure  3;  however,  the  overall 
values  of  the  scattering  rates  would  be  smaller  because  of  decrea.sed 
overlap  of  the  impurity  potential  with  the’  electron  wa\e  functions. 

4.  Random  Distribution  of  Impurities 

The  total  scattering  rate,  \V„  k.?'  <hie  to  a  macroscopic  number  of  ran¬ 
domly  distributed  impurities  can  be  found  from  the  one  impurity  scat¬ 
tering  rate  in  equation  (17)  by  summing  over  impurity  position  vectors 
Ro.  This  can  be  seen  by  considering  the  matrix  element  of  the  total 
potential,  which  is  a  sum  of  A’,  one-impurity  potentials,  each  cen¬ 
tered  about  a  position  vector  R^,  from  a  set  of  A’,  random  position^* 

'Hf  the  impurities  are  suKstitutional  then  these  vectors  coincide  with  crystal 
lattice  site.s. 
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vectors,  {Ro}.  The  matrix  element  for  scattering  from  N,  impurities  is 
then  a  sum  of  single  impurity  matrix  elements,  each  given  by  equation 
(11),  and  can  be  written  as 

(a,  exp(-iK  ■  Ro)/tpy(^J.  (20) 

{Ro} 

where  we  have  introduced  the  shorthand  notation  K  =  k'  —  k,  p  = 
{Q,k,(7}  and  p'  =  {o', k',^'}.  The  absolute  square  of  this  matrix  ele¬ 
ment  can  then  be  written  as 

|(o',k',q'lt/A/.|o,k,q}f  =  exp(-iK-(R,-R'J/lpy(z„)/t;p,(.;j. 

{Ro}  {r;,} 

For  a  macroscopic  number  of  impurities  at  random  positions,  the  phases 
in  the  sum  above  add  incoherently,  resulting  in  a  sum  that  is  close  to 
zero  except  when  R^,  =  R'^  or  K  =  0.  These  two  conditions  lead  to  two 
terms. 


Ko',k',^'lf/A.|a,k,^)p 


E  \^pA^o)\^  +  I  E  ^py(‘-o)r4.k'. 

{Ro}  {Ro} 


(21) 


VV 


o,k,7 


respectively,  where  the  sums  are  over  the  ^o-components  of  the  set  of 
random  impurity  positions.  The  second  term  is  negligibly  small  because 
the  random  phases  cancel  in  the  summation  over  Zo-  Substitiition  of  the 
first  term  into  Fermi’s  golden  rule  yields  a  total  scattering  rate  which 
is  a  sum  of  single  impurity  scattering  rates,  given  in  equation  (17). 

In  a  bulk  .semiconductor  the  single  impurity  scattering  rate  Wa.'k,q(~o) 
is  roughly  independent  of  position  Zg,  leading  to  a  total  scattering  rate 
proportional  to  impurity  density  7i,  =  N,f  L^.  However,  in  a  double¬ 
barrier  structure  the  function  can  have  a  strong  spatial  de¬ 

pendence  in  the  quantum  well  region  because  of  the  presence  of  the 
resonances.  For  on-resonance  (off-resonance)  electrons  the  wave  func¬ 
tion  amplitude  in  the  quantum  well  region  is  large  (small)  comf)arcd  to 
outside  this  region,  leading  to  a  large  (small)  contribution  to  the  total 
scattering  rate  from  impurities  in  the  quantum  well. 

The  sum  of  equation  (17)  over  impurity  positions  c„  can  be  ap¬ 
proximated  by  the  factor  n,  times  an  integral  over  the  crystal  volume. 
Integration  over  the  x„,yo  coordinates  contributes  a  factor  resulting 
in  a  total  scattering  rate  that  can  be  written  as 


(22) 
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where  the  integration  region  Ri  is  the  well  plus  barrier,  \zola\  <]  +  6, 
and  R2  is  the  region  outside  the  barriers,  1  +  (5  <  |^o/«|  <  A/Sa. 

The  first  term  in  equation  (22)  gives  the  contribution  to  the  total 
scattering  rate  from  the  quantum  well  (and  barriers)  region,  and  the 
second  term  gives  the  contribution  from  the  rest  of  the  crystal.  Since 
the  quantum  well  forms  a  negligible  volume  fraction  of  the  whole  crystal 
volume,  it  may  appear  at  first  sight  that  the  first  term  is  0{af  L)  times 
the  second  term,  and  hence  is  completely  negligible.  This  is  true  for 
electrons  that  are  in  states  away  from  the  resonance  regions.  However, 
electrons  with  wave  vectors  that  satisfy  the  resonance  conditions  (see 
sect.  1)  may  have  a  macroscopic  fraction  (i.e.,  not  0(nl L))  of  their 
normalization  in  the  quantum  well  region.*^  For  resonant  electrons, 
this  leads  to  a  situation  where  the  contribution  to  the  total  scattering 
rate  from  impurities  in  the  quantum  well  region,  given  by  the  first  term 
in  equation  (22),  is  comparable  to  the  contribution  from  the  rest  of  the 
crystal,  given  by  the  second  term.  This  is  essentially  a  finite  size  effect. 
As  expected,  for  sufficiently  large  crystals,  L  — >  00,  the  second  term 
always  becomes  larger.  We  find  that  for  a  sample  whose  total  thickness 
is  30  /i,  with  40-A  well  and  barriers,  and  an  impurity  potential  range 
1/A  =  10a,  the  first  term  is  comparable  to  the  second. 

5.  Conclusion 

Within  the  framework  of  a  simple  model,  we  explored  elastic  impurity 
scattering  for  electrons  in  resonant  states.  We  found  that  for  certain 
parameters,  such  as  screening  length  and  well  and  barrier  widths,  im¬ 
purities  in  the  quantum  well  region  can  contribute  significantly  and 
even  dominate  the  scattering  rate,  despite  the  fact  that  a  quantum 
well  occupies  a  negligible  volume  fraction  of  the  crystal.  We  have  used 
a  simple  model  of  isotropic  screening  which  is  an  oversimplification. 
Proper  treatment  of  screening  in  a  double-barrier  structure  is  in  itself 
a  difficult  problem.  We  have  also  ignored  effects  associated  with  space 
charge,  external  electric  field,  and  band  nonparabolicity. 
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Appendix  A.  Functions  in  Equation  (5) 


The  functions  in  equation  (5)  in  the  main  body  of  the  text  are 
defined  belo'v: 

=  ^exp(-u;(u))7;„(u) 

=  ^exp(u;(u))^^(u) 


/c(w)  =  ^exp(-f(l  +  (5)r/)|exp(ro(u)(5)  tj^iu)  +  exp{-iv{u)6)  '^a(^0} 

gdu)  -  -exp(i(l  +  6)u)|exp(?e(u)(5)  Vain)  +  exp(-u;(u)<5) 


7/e(u)  =  COs(u)  — 
r/o(u)  =  sin(ti)  + 
^e(u)  =  cos(w)  + 
=  sin(u)  - 


f3w{u) 


0w{u) 


sin(u) 


cos(u) 


u 


/3w{u) 


u 


sin(u) 

■  cos(u). 


0w{u) 

The  function  Fo,ic(u)  in  equation  (5),  correct  to  zeroth  order  in 
a/L,  can  be  written  as 

,2 


^a.k(w)  = 


u 


DM 


For  a  =  e,  the  function  Df{u),  in  the  denominator  is  given  by 


u'^) 

~2^w{ii)  «sin(2ri)sinh  {26iv{u)) . 

For  a  =  o,  the  function  Do{u)  can  be  obtained  from  D^{u)  by  the 
replacements  sin(u)  — >  cos(u),  cos(u)  — >  —  sin(u),  and  sin(2«)  — > 
—  sin(2u). 


sin  ^{u) )  (cosh  [28xr{u))  —  1 ) 


Z>e(u)  =  ^  (^7  +  (1  -  ^COS^(?i)  + 


/?(7  - 


Appendix  B.  Impurity  Matrix  Elements 


For  many  interesting  impurity  potentials  the  matrix  element  in 
equation  (8)  (see  main  body  of  text)  can  be  expressed  as  a  ma¬ 
trix  element  between  envelope  functions.  This  can  be  shown  by 
changing  the  integration  over  all  space  to  a  sum  over  all  lattice 
vectors  R  and  an  integral  over  a  single  primitive  cell  of  volume  v: 


{a',k',q\u\a,k,q)  =  k' ,,'(•*  + +  R)Wc(r  +  R).  {B  -  1] 

R  •''' 


The  envelope  functions  vary  slowly  over  a  primitive  cell,  so  they 
can  be  expanded  about  an  arbitrary  lattice  vector,  T/’o.k,(j(R-  +  r)  = 
+  ^’V’a,k,?(R-)  •  r  +  •  •  s  while  the  Bloch  functions  are 
periodic,  +  R)  =  Uc{r).  Using  this  in  equation  (B-1)  and 
keeping  only  the  first  term,  we  have 


(o',k',<7'|[/|o,k,9)  =  1:  /  d3rMr)l^t/(r  +  R)C-,kv(R)0.,k.,(R).  [B - 

R 


We  now  Fourier  expand  the  impurity  potential  in  plane  waves 
as  in  equation  (10)  of  the  text.  In  this  expansion  q  takes  the 
values  (Hj., Hy, nj)27r/U,  where  nj,,ny,nj,  are  integers.  Using  this 
expansion  in  equation  (B-2)  we  get 


{a',k\q'\U\a,k,q)  =  ^(q)  H exp(zq  •  (R  -  Ro))Vv.k',<,’(R)’/’a,k,<,(R)  /  d^r|ue(r)|^ exp(?q  •  r). 

q  R 


(B-3) 


If  wc  assume  that  f/(q)  %  0  for  |q|  >  </(-,  where  is  some  cutoff 
wave  vector  that  satisfies  q^  ■<C  l/|c|,  and  c  is  a  primitive  lattice 
vector,  we  can  cut  off  the  sum  on  q,  in  equation  (B-3),  at  )q|  = 

For  all  r  within  a  cell  we  have  q  •  r  <C  1,  so  we  can  expand  the 
exponential  in  equation  (B-3),  keeping  only  the  first  term,  and 
write 


I  d'V|i/c(r)|^exp(?q  ■  r)  =  v  +  0((7,|R1). 

J  V 

where  we  are  assuming  that  the  Bloch  functions  are  normalized  to 
V  over  the  primitive  cell.  Using  these  approximations  in  equation 
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(B-3),  changing  the  sum  on  R  to  an  integral  ov'cr  the  crystal 


volume  f>, 


(R-5) 


we  can  write  the  matrix  element  in  equation  (B-3)  approximately 
as  given  in  equation  (9)  of  the  text.  Note  that  a  delta-function 
impurity  potential  does  not  allow  this  simplification. 
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Appendix  C.  Definition  of  the  Function 

In  this  appendix  we  give  the  explicit  form  of  the  quantity  Ia'.k',q'-,a.k.g(~o)- 
introduced  in  equation  (15)  of  the  text.  When  the  impurity  is  in¬ 
side  the  quantum  well  region,  0  <  Zo  <  a, 


'a'.k'. 


9  ,a,j 


,(^o)  =  T 


4[D^.{q'a)D,{qa)Y/^ 


aqq  kc+o.'  [exp(-sro) 


exp(-s  -f  i[q  -f  q'))  -  1 
-s  i[q  +  q') 


+  (-1)"  exp(-szo) 
+  (-l)“+“'exp(s-^,) 

+  (-l)"exp(,srj 


exp(-s  +  i{q  -  q'))  ~  1 
-s  +  i{q  -  q') 

expj-s  +  tjq  +  q'))  -  exp[{-s  +  ijq  +  q')): 
—s  -f  i{q  +  q') 

exp(-s  +  iiq  ~  q'))  -  exp[(-s  +  i{q  -  q'))z^ 


-s  -f  i{q  -  q') 

-  (-l)^exp(-s2o) 


s-^i{q  +  q') 

1  -  exp[is  +  i{q  -  q'))z^ 


s  +  iiq  -  q') 


+  ^[exp(s(:;o  -  1))  +  ( exp(-s(2o  +  1))]  [v(g'a)  5?a(ga 


cxp[(— s  +  A  +  A'*)5a]  —  1 
'  -,s  +  A  +  A'* 


+  I'Aq  a)  Vaiqa) 


exp[(— s  +  A  —  A'*)A(7]  —  1 
-s  -f  A  -  A'* 


,  ,pxp[(-.s  -  A -f  A'*)6a]  -  1  ^  exp[( -s  -  A  -  A'*  )(5n] 

+  qAq  (i)iaiq(i) - ,  ,  - +  L-iq 


-.s  -  A  4-  A'* 


-.s  -  A  -  A'* 


+  fAq'(^)  foiqn) 


exp(— S(5a) 


^  /.(?«) 

-f(q-q') 


cxp(  — .sAn) 
s  -  i(q  +  q') 


,  \  ^  N  Pxp(-s6n)  ,  ,  exp(-s(5a) 

+  fc‘{q  o-)9c,{qa)—-— — —7+  ga'(q(i)9c.{qn]——, - rrr 

s  +  i{q+q')  s  +  i{q-q')i^ 


iC-l) 


2.3 


When  the  impurity  is  inside  the  barrier  region,  a  <  r,,  <  ( 1  + 


■^o',k',i7';a,k,g(^o)  —  . 


4  s^o)R.+a’ 


exp(-.s  +  i{<j  +  q'))  -  1 

-s  +  i(q  +  q’) 


+  (-l) 


cxp{-s  +  t(q  -  q'))  -  1 


-.s  +  t{q  -  q') 


+  (-!)"+"' exp(-s.~„) 

+  (-l) 


exp(s  +  ii(?  +  g'))  -  1 
6  +  i(7  +  q') 

exp(s  +  i{q  -  q'))  -  1 


s  +  iiq  -  q') 


1  • 


«  '■  — xS  +  A  +  A  ' 


+  ''L'W(i)Liqo.) 


exp[(-s  +  A  -  A'')a<S]  -  1 
— +  A  —  A  * 

cxp[(-,s  -  A  +  A’')rt^]  -  1 
-  A  +  A'* 


+  C’W(i)L{q(i) 


exp[(-,s  -  A  -  A'*)n^]  -  1 


-5  -  A  -  A'* 


+  -  <’xp( 1))  fi^.{q'a)T)„{qa) 


exp[(As  +  A  +  A'*)(c^  -  1 )]  -  1 
.s  +  A  +  A'* 


+  Oiq'^)vciqfi) 


exp[(.s  +  A  -  A'*)(::,  -  1 )]  -  1 
s  +  A  -  A'* 


+  q'a'{q'a)(.4qa) 


^exp[(s  -  A  4-  A'-)(Co  -  1 )]  -  1 


•s  -  A  +  A'* 


+  Oiq'fOLiqa) 


exp[(.s  -  A  -  A'* )( Co  -  1)]  -  1 
s  -  A  -  A'* 


24 


+  -  exp(.s(-„  - 
a 


n'A(i'(i)ic{qa) 


exp[(-5  +  A  +  A'*)fl^]  -  exp[(-.‘;  +  A  +  A'*)(c„  -  1  )| 
-s  +  A  +  A'* 


exp[(  -s  +  A  -  A'’)a<5]  -  exp[(  -.s  +  A  -  A'*  )(co  -  1 )] 
— 5  -p  A  —  A'* 


+  '/.V  (<■/«)  d-.  (<■/«) 


exp[(-.s  -  A  +  A'')a^]  -  exp[(-.‘;  -  A  +  A’')(r^  -  1 )] 
-.s  -  A  +  A'* 


oxp[(-.s  -  A  -  A'*)a6]  -cxp[(-.s  -  A  -  A'*)(c^  -  1 )] 


-5  -  A  -  A'* 


-  -  Pxp(--"«(1  +  ^))(expt^5::o)  +  (-ir‘^'''('xp(-.^cj)  [-^ 


A'jq'd)  Iri go) 

+  ^ig  -  g') 


9a'ig'»]  loigf^)  ly’jg'fi)  9c,{gfj) 

-■■<  +  iig  +  g')  -s-)(q  +  q')  -s-i{q-q')iy 

In  equations  (C-l)  and  ((’-2)  vve  have  introduced  the  functions 
gdgo)  =  qari-,(qa). 


L{g(‘)  = 

faigo)  =  '2g(i  nxp[/<ifrt((5  +  1  )]/„(?«), 

(ja{qa)  =  2f/oexp[-?9fl(^  +  1  )].(7,^((7a). 

A  =  —xr{qn). 
a 

where  r/„(<7n),  ^0(70),  faign)^  7n(7«)-  and  ii'{u)  are  defined  in  ap¬ 
pendix  A  and  equation  (6),  respectively.  When  these  functions 
appear  with  primed  sub.scripts,  the  final  state  value  of  the  trans¬ 
verse  wave  vector  (=  k')  should  be  used  in  their  definition.  The 
parameter  A',  for  the  final  state,  is  defined  by  A'  —  u'{q'a)/a. 
where  k'  is  used  in  equation  (7)  in  place  of  k.  We  have  also  in¬ 
troduced  the  operator  /?„,  which  operates  on  a  complex  number 
r  according  to  the  definition  R„  :  =  c  +  (  —  1)";*,  where  a  is  an 
integer.  In  equations  (C-l)  and  (C-2)  we  use  the  convention  that 
o  =  0  for  an  even  state  and  rv  =  1  for  an  odd  state. 
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